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The transport of quantum electrons through hierarchical lattices is of interest because such lat- 
tices have some properties of both regular lattices and random systems. We calculate the electron 
transmission as a function of energy in the tight binding approximation for two related Hanoi net- 
works. HN3 is a Hanoi network with every site having three bonds. HN5 has additional bonds added 
to HN3 to make the average number of bonds per site equal to five. We present a renormalization 
group approach to solve the matrix equation involved in this quantum transport calculation. We 
observe band gaps in HN3, while no such band gaps are observed in linear networks or in HN5. 

PACS numbers: 64.60.-i, 02.70.-c, 02.70. Ns 

Keywords: Hanoi network, renormalization group, matrix RG, quantum transmission 



I. INTRODUCTION 

Understanding and controlling the transport of elec- 
trons is central to the operation of all electrical and elec- 
tronic devices pQ . In many cases of interest in nanoma- 
tcrials, the electron transport is coherent, and therefore 
must be analyzed using the Schrodinger equation. Inter- 
ference effects can then lead to a metal-insulator transi- 
tion and Anderson localization [2]. Even more than fifty 
years after Anderson's publication of his celebrated pa- 
per, the effect remains an active area of study [3j-[9] . The 
main goal is to calculate the transmission probability as 
a function of energy, E, for an incoming electron, i.e. the 
probability that an electron that comes in from x = — oo 
can be observed at x = +oo. 

The starting point to quantum calculations of (spin- 
less) electrons through a material is often a tight-binding 
model [lfj . In such a model each node can be consid- 
ered an atom, the on-site energy at a node is associated 
with a potential energy at the site, and there is a hopping 
term which comes from discretization of the kinetic en- 
ergy term of the Schrodinger equation |10j . The electron 
transport calculation via the Schrodinger equation thus 
reduces to the solution of an infinite matrix equation. 
The solution of the matrix equation is often accomplished 
using a Green's function method O H~0hT5] . In this paper 
we instead use an ansatz approach introduced by Daboul, 
Chang, and Aharony |16j , which is simpler to describe at 
an undergraduate level [T7] . The ansatz reduces the size 
of the matrix equation to that of the number of tight - 
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binding cites in the scattering volume, plus one for both 
the incoming lead and the outgoing lead. This approach 
has been used by other authors [18] . We find that the 
ansatz approach is particularly well suited to our calcu- 
lations of transport through hierarchical structures. To 
perform the calculation we construct a decimation Renor- 
malization Group (RG) procedure to reduce the number 
of sites of the hierarchical structure. Our RG is related 
to that utilized by others for tight-binding models [T51 - 
[22] • However, our RG has been explicitly constructed for 
the calculation of the transmission probabilities. We find 
that our RG procedure greatly simplifies the calculation, 
albeit only for certain select networks that have hierar- 
chical structures. Although our RG does not significantly 
simplify the calculation of associated wavefunctions, we 
nevertheless give a recipe for the RG calculation of the 
wavefunctions. 

The specific models we solve here for quantum trans- 
port are motivated by four considerations. One is to un- 
derstand how hierarchical models, in particular the Hanoi 
networks [23) . affect quantum transport. The second is 
that such hierarchical models provide an intermediate be- 
tween regular lattices and ones that have a small-world 
property [24], and would be of interest to understand 
quantum transport of nanomaterials that have a small- 
world property [5] !25fE5] . The third is that often phase 
transitions such as the metal-insulator transition or a fer- 
romagnetic transition have some universal quantities that 
depend only on the dimension. Hierarchical models then 
sometimes provide insight into how these universal quan- 
tities behave as a function of the dimension [201 I21)rt3"2] . 
Lastly, experimental realizations of hierarchical materials 
may possess novel physical properties [3"5rl3"T] . 

The hierarchical models we study are the Hanoi net- 
works HN3 and HN5. These networks have particularly 
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interesting properties. First, they are both planar, and 
consequently could be experimentally constructed on a 
surface. Second, both networks have typical 'paths' (de- 
fined precisely in Sec. [TlJ) that grow more slowly with 
system size than do paths in regular lattices. Finally, 
Anderson localization is associated with randomness in 
the system, while randomness in the Hanoi networks de- 
pends on the scale. For example, for HN3 locally every 
site has three bonds connecting it with other sites, choos- 
ing sites from HN3 at random the connections to other 
sites seem random, but these connections are actually 
from a hierarchical arrangement and hence at the larger 
scale there is regularity to the network. Therefore study- 
ing electron transport and Anderson localization in these 
lattices is of interest. 

In Section|ll]we provide a brief description of the Hanoi 
networks. For transport properties, these networks can 
be connected to the leads in many ways, but we choose 
to present results only for symmetric linear and symmet- 



ric ring lead attachments. In Section [TXT] we develop the 
RG equations for calculating the transmission through 
these networks, with details of the RG presented in Ap- 
pendix [XJ In Section IV we analyze these RG equations 
for these networks. This involves iterating the RG until 
the system is comprised only of a few lattice points, and 
these small lattice solutions are presented in Appendix |B| 
SectionlVlcontains our conclusions and a discussion of our 
results. The appendices contain the basic matrix algebra 
used to develop the RG equations. 



II. NETWORK STRUCTURE 

Each of the networks considered in this paper possess 
a simple geometric backbone, a one-dimensional line of 
N = 2 k sites formed into a ring as depicted in Fig. [I] 
for HN3 and k = 5. Alternatively, we can connect the 
network to the incoming and outgoing leads in a linear 
arrangement with 2 k + l sites, as depicted in Fig.[2]and|3] 
Each site is at least connected to its nearest neighbor left 
and right on the backbone. For consistency, we call the 
ordinary one-dimensional ring HN2 (for Hanoi Network 
of degree 2). For example, HN2 is the linear lattice in 
Fig. [3] formed by only the black bonds. 

To generate the small-world hierarchy in these net- 
works, consider parameterizing any number n < N (ex- 
cept for 0) uniquely in terms of two other integers 



i > 1 and 1 < j < 2 k 



= T- 1 (2j-l). 



(1) 



Here, i denotes the level in the hierarchy whereas j labels 
consecutive sites within each hierarchy. To generate the 
network HN3, we connect each site n — 2 l ~ 1 (4j — 3) 
also with a long-distance neighbor n' = 2*~ 1 (4j — 1) for 
1 < j < 2 k ~ 1 - 1 . (In the ring, if an index n equals or 
exceeds the system size N, we assume that the site n 
mod N is implied.) 




FIG. 1: Depiction of the 3-regular network HN3 with a one- 
dimensional periodic backbone forming a ring, here with k = 
5. The top and bottom sites n — and n = 2 fc_1 require 
special treatment and are connected to external leads. With 
these connections, the entire network becomes 3-regular. Note 
that the graph is planar. 

For the linear arrangement, the sites zero and 2 k + l are 
connected to the input/output leads, and site 2 fc_1 will 
not be connected to any other site or to the input/output 
leads [see, for instance, HN5 in Fig. [2], For the ring 
arrangement, the sites with number zero and 2 k ~ 1 are 
connected to the input /output leads [see Fig. [I], and site 
zero is connected to site 2 k — 1 to form the ring. 

Previously [53], it was found that the average "chemi- 
cal path" between sites on HN3 scales as 



iHN3 



(2) 



with the distance I along the backbone. In some ways, 
this property is reminiscent of a square-lattice consisting 
of N lattice sites with diagonal ~ y/N. 

While the preceding networks are of a fixed, finite de- 
gree, we can extend HN3 in the following manner to ob- 
tain a new planar network of average degree 5, hence 
called HN5, at the price of a distribution in the degrees 
that is exponentially falling. In addition to the bonds in 
HN3, in HN5 we also connect all even sites to both of its 
nearest neighboring sites within the same level of the hi- 
erarchy i > 1 in Eq. 0. The resulting network remains 
planar but now sites have a hierarchy-dependent degree, 
as shown in Fig. [2j To obtain the average degree, we ob- 
serve that 1/2 of all sites have degree 3, 1/4 have degree 
5, 1/8 have degree 7, and so on, leading to an exponen- 
tially falling degree distribution of V {a = 2i + 1} oc 2~ l 
for i = 1, 2, 3, • • • . Then, the total number of bonds L 
in the (linear) system of size N = 2 k + 1 is 



fe-i 



2L = 2{2k - 1)+ V (2i + 1) 2 k ~ l = 5 x 2 k - 4, (3) 
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FIG. 2: (Color online) Depiction of the planar network HN5. 
Black lines demark the original HN3 structure, the green- 
shaded lines are added to make HN5. Note that sites on the ^ 
lowest level of the hierarchy have degree 3, then degree 5, 7, 
• • • , making up a fraction of 1/2, 1/4, 1/8, • • • , of all sites, 
thereby making for an average degree 5 of this network. 




FIG. 3: (Color online) Scattering a quantum electron off a 
linear version of HN3, here drawn as a branched Koch curve. 
The incoming electron on the left gets scattered into an out- 
going transmitted portion (right) and reflected portion (left) 
on the attached external leads (blue-shaded). In this form of 
HN3, the one-dimensional backbone is marked by black links 
while the small-world links are shaded in red. Note, for in- 
stance, that the shortest end-to-end path here is the baseline 
of the Koch curve. 



thus, the average degree is 



2L 
~N 



■5. 



(4) 



In HN5, the end-to-end distance is trivially 1, see 
Fig. [2} Therefore, we define as the diameter the largest of 
the shortest paths possible between any two sites, which 
are typically odd-index sites furthest away from long- 
distance bonds. For the N — 32 site network depicted 
in Fig. [2j for instance, that diameter is 5 as measured 
between site 3 and 19 (0 is the left-most site), although 
there are many other such pairs. It is easy to show re- 
cursively that this diameter grows strictly as 



K k 



FIG. 4: (Color online) Decimation of a block in the RG. The 
two sites with on-site energy «i with a connecting bond of 
strength T\ are decimated. 



III. MATRIX RG FOR HANOI NETWORKS 



At each RG step, we decimate all the odd sites. Take 
site to be at level i = k + 1 in a linear geometry and 
at level i — k in a ring geometry. As the odd sites each 
have only one small-world-type bond, we can divide the 
network into blocks containing 5 sites and decimate the 
pairs of two odd sites block by block. Let us start with 
the first block which contain sites 0,1,2,3 and 4. This dec- 
imation process for a linear geometry is shown in Fig. [4] 
Thus, here (see Appendix [A| we have 
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,HN5 



2 [k/2\ + 1 - log 2 N 



(5) 



After decimation, 



with [x\ the integer portion of x. We have checked nu- 
merically that the average shortest path between any two 
sites also increases logarithmically with system size N. 
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r k, 1 t 



A BD^BT. 



(8) 
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After simplification, we find that 



K 2 



a; 



«' 2 



«3 



2r 2 

ki + n 

r 2 Ki 



= Ai- 



A 2 + 



«1 + Ti 



(9) 
(10) 

(11) 
(12) 
(13) 



Here, the primed and unprimed quantities represent the 
1-st and the 0-th level respectively in the RG recursion. 
Notice that the decimation of sites connected to site 4 
[the right-most site in Fig. |4^a)] is not complete yet and 
therefore its on-site energy will be modified further when 
we decimate the odd sites of the next block which contain 
sites 4,5,6,7 and 8. 

After the decimation of the next block we get 



«2 



t 2 ki 



«3 



2t,?ki 



(14) 



Continuing the decimation block by block in this way, we 
find that 



K 



i+1 



2t 2 ki 



A' 



n+i Vi > 1 
Ai+i Vi > 2. 



for Vi G {2, . . . , k} (15) 

(16) 
(17) 



At first it appears that there are a lot of RG variables 
to worry about. However most of these RG variables 
are interdependent. It can be deduced from Eqs. ( 9|15 
16|17| that 



r ( m ) 
l 

Am) 



Tm+l 



= A 



m+2) 



(18) 
(19) 



and that the on-site energy parameter of the even sites 
is related to those of the odd sites as 



Am) 



+ K m+l - K m +i - 2(A 1 m) - A m+ i) for Vi G {2, . . . ,k — m} 



and that specifically for a linear geometry, the on-site energy parameter of the end sites 



K fe+ i + (4 m) - n m+2 )/2 for Vme{0,...,*-2}. 



(20) 



(21) 



Thus we are left with just three independent RG variables which are k\ , Tq and A^ governed by the RG equations 



Am+i) _ (m) 



K-m+2 — Hm+1 — 2(A^ ^ — A m -|_l) — 



_(m+l) _ ,(m) 



'0 



= Ai" 



2[r (m) ] 2 

K l T Trn+1 

(m) , 
K l T Trn+1 



r ( m )l2 

>(m+l) _ x . L r J T m + 1 

^1 — ^m+2 T 



„( m )l 2 _ _2 
"a J T m+1 



Vm e {0, . . . , k - 2} 
Vm e {0,. .. ,k - 1} 
Vrne {0,...,fc- 1}. 



(22) 
(23) 
(24) 



r 



IV. ANALYSIS OF THE RG EQUATIONS 

A. One-dimensional Lattice (HN2) 

It will prove helpful to demonstrate the general set 
of recursions [Eqs. (22 23 24 1] by way of the one- 
dimensional (d = 1) ring of N — 2 k sites. For consis- 
tency we call this the HN2 network, or Hanoi network 
with 2 bonds per site (with additional bonds for the in- 



put/output leads in the ring geometry). We can em- 
ploy the recursions to explore the transmission through 
a d = 1 ring of N = 2 k sites. The energy scale chosen 
throughout is such that the (uniform) transmissivity for 
each bond has a unit weight. With that, we obtain the 
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initial conditions 



"i — 

.(0) _ 

— 

JO) = 

i 

v(0) _ 



E, (i > 1) 
-1, 

0, (i>l), 



(25) 



0, 



Eqs. ([22] [23] [24| simplify to 



(i > 1) 
2t 2 



(26) 



'rn+l — 



Am) 



and t, = A, = for all i > 1 and 



where K m 

— ( m ) 
'm — 

by denning 

obtained by dividing the 2 nd by the 3 rd line in Eqs. (261 
Formally, the solution is 



These nonlinear recursions are easily solved 
j/t ot for which s m+i 



2 cos 



2r(°) 



2T 2 



2t(°) 



(27) 

where T n (x) refers to the n-th Chebyshev polynomial of 
the first kind [35]. Inserting into Eqs. (26) and applying 



the initial conditions in Eqs. (25), generates the results 

-l 



_(m) _ 



Am) 



m— 1 1 

■n 



(28) 



i=0 



where the last equality emerges under reordering factors 
in the products. 

Eq. ( B9 ) in Appendix [B] shows that the transmission 
amplitude t is directly proportional to . Clearly, if 
there is no transmission on any bond, i.e. r( fe ) — 0, for 
a given input energy E, there can be no transmission 
through the network itself, no matter what happens on 
the sites. But instead of plotting r^ k \ it will prove more 
instructive to plot K^ k \ It is easy to see from Eqs. (28) 
that «( fc ) varies rapidly whenever does, but that k ( 1 
varies smoothly whenever vanishes. In the following, 
we will see that this behavior remains true for HN3 and 
HN5, in which case the variation of kP*' with = E 
provides more information beyond the mere vanishing of 

Uk) 



We have evolved the RG-recursion in ( 26 1 for the ini- 
tial conditions in (25) and plotted K t - k = TU > as a func- 
tion of k(°) = E in Fig. [H] Even at that system size, 
N = 2 k = 1024, delocalized states completely cover the 
domain -2 < E < 2. 



Case HN3 



Again, we can employ the RG Eqs. (22 23 24) for 
transmission through HN3 consisting of a ring of N = 2 fc 




FIG. 5: (Color online) Plot of k (10) in Eq. (Q as a function 
of its initial condition = E. Even at this small order, the 
function varies extremely rapidly, such that its (green-shaded) 
line completely covers the shown domain. Thus, for —2 < 
E < 2, K (m) for any sufficiently large m is a random function. 
Correspondingly, the transmission spectrum is dense, with full 
transmission close to any input energy E such that particles 
do not localize. 



sites, as in Fig. [T] Since all initial diagonal entries are 
identical, the hierarchy for the collapses and we re- 
tain only two nontrivial relations, one for k± and one 
for all other ft, = K2 for all i > 2. Here, all n are 
non-zero, encompassing the backbone links (i = 0) and 
all levels of long-range links (i > 1). But it remains 
t,; = —1 for i > 1 at any step m of the RG, in par- 
ticular, = — 1 throughout; only the backbone tq 
renormalizes non-trivially. Although all links of type Aj 
are initially absent in this network, the details of the 
RG calculation shows that under renormalization terms 
of type Ai emerge while those for Aj for i > 2 remain 
zero at any step. Thus, we obtain far more elaborate RG 
recursion equations compared to those of HN2. Abbrevi 



ating K n 



- J m ) 



- J™) 



and A m = A^ m) , Eqs. ([22 



23 24 ) and their initial conditions reduce to 



Am+1 



2A r 



(kq = E), 
(to = -1), 
(A o = 0). 



(29) 



We have evolved the RG-recursion in ( |29| and plotted 
^^=200 as a function of kq = E in Fig. [6] Even at 
that enormous (and definitely asymptotic) system size, 
N = 2 200 pa 10 70 , domains of localized states remain 
asymptotically inside the physically relevant domain of 
— 2 < E < 2. For comparison, the radius of the visible 
universe is only about 10 fm. 

In the next subsection, we will explore the asymptotic 
properties of these recursions for large m. We will find 
domains in E of stationary solutions, which are particular 
for HN3, and show that the special points where this 
analysis fails correspond to transitions between localized 
and delocalized behavior. 
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FIG. 6: (Color online) Plot of K200 in Eq. ( |29[ ) as a function 
of its initial condition kq = E. Bands of localized and delocal- 
ized states intermix. Correspondingly, there are localization- 
delocalization transitions already before the addition of any 
additional randomness in HN3, merely as a function of the 
input energy E. On t he h orizontal axis, we have marked the 
solutions E { P of Eq. Q for s = 1 (red dots), s = 2 (blue 
dots), and a = 6 (small black dots). The accumulation of the 
latter demonstrates (even for such a small value of s) that 
the band gaps are associated with the absence of such solu- 
tions. While the solutions for s = 1 happen to be interior to 
the bands, some of those for s = 2 appear to mark the band 



edges, in particular the one at E = E. 



-0.637875. 



where C > is another undetermined constant that de- 
pends on E. Using the recursion for n m in Eq. (29) to 
next-to-leading order yields 

K m+1 ~ K m 2A m , (33) 

~ K m + 2(^-l)exp{-CV2^ T }, (34) 

which, when summed from a m > too to 00 results in 

n m ~ Koo - 2 (/& - 1) exp |-CV2'»+ 1 } , (35) 

where we have kept only the first term in the sum, as the 
summand is exponentially decaying. 

The consequences of this analysis are quite dramatic. 
If such a steady-state solution is reached, both transmis- 
sion rates r and A vanish, only leaving a finite on-site 
energy Kqo- Hence, there can not be any transmission 
through the network when such a state is reached, and 
gaps emerge in the transmission spectrum. Since the sys- 
tem size is given by N = 2 m , this result implies that finite 

size corrections scale with exp /— C\/2N X, i.e. finite-size 

corrections decay rapidly with a stretched exponential. 
Note that there are no steady-state solutions that cross 
= ±1 (dashed lines in Fig. |6|, where the correction 
in Eq. (35) would break down. Instead, the approach of 



Koo — > ±1 frequently appears to be associated with the 
emergence of a band edge between localized and delocal- 
ized states. 



Analysis of the Steady State: 



We can analyze the absorbing steady state, which is 
the unique feature of HN3 (in contrast to HN2 and HN5) 
leading to band gaps, as follows. Numerical trails show 
that the RG recursions in Eq. ( 29 1 reach a steady state 



for certain initial conditions E when for all m larger than 
some mo it is 



1 > A™ > 



K m 1 



(30) 



The leading contribution for K m in Eq. ( 29 ) then suggests 

K m+1 ~ K m ~ K c»! (31) 

which is a constant that is difficult to derive from the 
initial conditions, kq — E, unfortunately. 

From the recursion for r m in Eq. ( 29 ) , we further ob- 
tain 



T m +1 



A.: 



'm-1 
2 - 1' 



(32) 



This 2nd order difference equation has two solutions, of 
which we discard the oscillatory one, to get for large to > 
m : 



1) exp{-CV2™j , 
(4,-1) exp{-CV2^} 



2. Band-Edge Analysis: 

The non-trivial band structure warrants some further 
investigation. In particular, we can associate such band 
edges with initial conditions E = E^ for which there 
exists a m — s such that 



« s = -1, 



(36) 



a singular point in the recursion for A s +i. [Interestingly, 
any singularity at K m = +1 appears to be benign in 
that it does not affect the continuity m /?™ cis Sl function 
of E for to — > 00; it afflicts each quantity in Eqs. (29) 
simultaneously, leading to a divergence in K m +ii Tm+ij 
and A m _|_i just so that n m+ 2 ~ K m-i, T m+2 ~ T m-i> 
A m +2 ~ A m _i.] In Fi g, pi we have also marked the real 
solutions E^ of Eq. ( |36| for s — 1, 2, and 6. Clearly, 
those solutions strongly correlate with the bands, and 
there appear to be none within the gaps (although we 
have not been able to prove this conjecture). But while 
those solutions for 8 = 1, E^' 2 ^ = ±a/3 7 are located well 
within some band (as are those for s = 6), the four real 
solutions of Eq. ( 36 1 for s = 2 satisfy the quartic equation 



= 1 - 2£ - 6£r 



E 4 



(37) 



and appear to be all associated with some more or less 
significant band edge, see blue dots in Fig. [6] We can 
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speculate that there is a whole hierarchy of transitions, 
each associated with one of the solutions E rn ! , which 
may become dense on certain intervals. While we don't 
know what determines those intervals precisely, we can 
analyze the behavior in the neighborhood of Eq. ( 36 1 . 



We observed that the recursions in Eq. ( 29 ) possess sta- 



ble steady-state solutions for large m characterized by 
r m ~ A m —¥ 0, i.e. vanishing bond-strength between 
input and output. These solutions prevail in the ob- 
served band gaps, which accordingly correspond to local- 
ized states. It seems that the reason for the persistence 
of gaps derives from that stability: band gaps emerge 
whenever the steady state is reached before Eq. ( 36 1 can 



be satisfied. For instance, in the case of HN5 below, any 
putative steady-state solution proves unstable for suffi- 
ciently large m such that any band gaps are transitory 
only, see Fig. [9| 



For the analysis of the recursion Eq. ( 29 ) , we assume 
that for some m = s, we reach 



-1 



(e<l), 



(38) 



where e = e(E) may be of either sign, depend ing on 



AE = E - E. 



Generically, e oc AE, see Sec. 



below. Assuming that r s , A s <C 1/e leads to 



IVB3 



T s + 1 
A s +1 



-1-2A, 



0(e), 



2c 



4 



0(e), 
"0(e), 



(39) 



which leaves only A s +i singular. After one more recursion 
step, we get instead 



Ks+2 
T s +2 



.11 

e 



2 + 6A S + 2A 2 - 2t 2 - 5A s t 2 



-2 - 2A, 



2e 4 
A s+2 ~ 0(1). 



r 2 + 3A s r, 



2 

s's 



4A e - 2r 2 



0(e), 



0(e), 
(40) 



At this point, Eq. (29) decouple to leading order, as 
A s +2+ t +i ^ +2+ Jl? s+2+l 1/4 remains of order 



0(1) while both n s+ 2+i and r s+ 2+i are of order 0(1 /e), 
and we get for some i > 



Ks+2+i + l 
Ts+2+i+l 



Ks+2+i — 

' s+2+i 
K s +2+i 



2t 2 

ZT s+2+i 

; 

K-s+2+i 



(41) 



These are exactly the same recursions we obtained in 
Eq. J26| for HN2, with the solution in Eq. p7): 



where from Eq. ( 40 ) we have 

Ks+2 



2t. 



s+2 



A = 



-1+A^r, 



2 + 6A S - 3t 2 - 8A s t s 2 
2 + 2X - t? 



(43) 
(44) 



With that inserted into Eq. (42), we can deduce 



i'.s+2 + f 



2r. 



s+2+i 



Toi 1 - A- 



1 - A-jF'v (1) , 



1 - 2 ll A- 



(45) 



since T' n [x) — nU n -\(x) and U n -\(1) = n, referring to 
the Chebyshev polynomial of the 2nd kind, U n [x) |38j . 
With the exponential growth in i of the correction ampli- 
tude in the asymptotic expansion in Eq. ( |45[ ) , the expan- 
sion breaks down at some i ~ Iq such that the correction 
itself becomes of 0(1), i.e. 



'o 



1 , ( T s 

2 l0g2 \W\ 



(46) 



For i > i a , according to the first line of Eq. ( 45 ) the ratio 



n s +2+i/T s +2+i either rises or falls exponentially, depend- 
ing on whether Ae < or Ae > 0, respectively. In the 
latter becomes less relevant and the bonds r and 

A determine the future evolution in m, leading again to 
the chaotic behavior in n m observed within the bands 
in Fig. [6] On the other hand, if Ae > 0, the on-site 
energies K dominate exponentially over the couplings t 
and A, evolving towards an absorbing steady state on the 
band-gap side of the transition. 



3. Scaling Relation for Koo (E) : 



As mentioned in Sec. IV B 1 we can not generally pre- 



dict the dependence of the asymptotic behavior on the 
initial condition E. But we can use the analysis of that 
section at least to determine the behavior of k oq (E) on 
the approach to those band edges where it diverges. 

It is easy to show, using just the singular terms for k 
and r in Eq. (40) inserted into the recursions in Eq. (41 ), 



that both simultaneously decay exponentially with i 
least while i < io from Eq. 



at 



(46), 



K s +2+i 
T s +2+i 



2 l e 
2*+ 1 e 



0(1), 
-0(1). 



(47) 



Ks+2-H 

2^ 



s+2+i 



Hs+2 



It. 



s+2 



(42) 



At i ~ i there is a cross-over beyond which for all i > 
io it is T s+ 2+; and a saturated value of n s +2+i 



C. Interpolation between HN2 and HN3 




yj E -E' 



0.02 0.04 0.06 0.08 0.10 



FIG. 7: (Color online) Plot of 1/noo as a function of 

r(2) _ 



« for E -> Ei l) 
0.637875, marked as the 2nd blue dot from 



E - E\ ' for E -> EV to test Eq. \49\. Here, E\ 



the left in Fig. [6] 



Koo is reached. We can obtain the dominant asymptotic 
behavior at the cross-over from 




(48) 



We can further establish a (generic) relation between e 
and 



AE ~ E-E® 



by extending the discussion of Eq. (36) in Sec. IV B 2 
We set 



~ « s (^) + <(^ i3 ) 



AS, 



AS, 



since k s — > — 1 is a regular limit for AE 1 — > 0, as Eq. (37), 
for example, suggests. Hence, 



:(A£) ~ 



AS, 



and we conclude 



14" 



sgn k' E, 



)AE) 



with AE — > 0. In Fig. [7 
totic relation 1/kqo ~ y 



(49) 

AS| 5 

we have tested the asymp- 



AE\ in the band gap near 



E 



(2) 



in Fig. [6j 



-0.637875 . . ., a solution of Eq. (37) marked blue 



It proves fruitful to consider an interpolation between 
the case of HN2 in Sec. HVAl and HN3 in Sec. HVBl in 
terms of a one-parameter family of models. To wit, we 
can accomplish such an interpolation by weighting the 
transmission along the small-world links (see red-shaded 
links in Fig. [3]) by a factor of y relative to that of the 
backbone links (see black links in Fig. [3]). Clearly, more 
generally hierarchy and/or distance-dependent weights 
could be introduced as well. For y = 0, small-world links 
are non-existent, and we have the linear lattice HN2. Al- 
though still mostly delocalized, the states of the systems 
immediately change behavior when y > 0, and we find lo- 
calized states which expand their domain towards y = 1, 
corresponding to HN3, and continue to do so until at 
about y = 4 no transmission is possible any longer: The 
more we weight small- world links here, which classically 
would expedite transport 39 , the less quantum transport 
is possible! In the next section we will see that even more 
small- world links, as in HN5, can lead to more transmis- 
sion again. Hence, the detailed structure of the links 
matter. 

To explore this y-family of models, we have to gener- 
alize Eq. (29) appropriately: 



— 2A„ 



•2t;, 



Kn+1 — 



V T rr 



(k =E), 

(to = -1), 
(A o = 0), 



(50) 



since for all non-backbone links in Eqs. ( [22| [23| [24] ) it is 



.(0) 



(0) 



-y at every RG-step. Note 



that these equations reduce to Eqs. (26) for y = (with 
all \i ee 0) and to Eq. ^ for y = 1. 

In Fig. [8] we map out the state of Eq. (50) after the 
1000 th iteration based on whether a steady state has been 
reached or not, depending on the incoming energy E and 
the relative weight y. For any y > 0, the ability to trans- 
mit has a strong chaotic dependence on these parameters, 
and ceases completely for y > 4. Even within domains of 
apparent transmission there are often sub-domains where 
no transmission is possible, and it is not clear whether 
true conduction bands exist. Since in this model the 
long-range links are not connected to each other except 
through the backbone, one may speculate that even at 
high weight these links merely lead to localized reso- 
nances that interfere with transport along the backbone 
instead of conveying it. (A similar confinement effect was 
observed for the RG applied to random walks on HN3 in 
Ref. [40].) 

It is straightforward to generalize the discussion for 
HN3 in Sec. |IVB| to this model. In particular, for the 
band-edge analysis we have to generalize Eq. ( 36 1 to read 



-y- 



(51) 



For s = 0, 1, 2, and 3, the numerical solutions E^\y) are 
also plotted as lines in Fig. [8] The result underlines the 
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contention made before for HN3 that the transitions be- 
tween transmission and localization are closely associated 
with these singular points of the Eq. (f50|). For instance, 



the big blue-shaded dots in Fig. [6] correspond here to the 
intersection of the simple-dashed line for s = 2 with the 
dotted horizontal line along y = 1 (i. e. HN3). Dominant 
features emerge, such as the line Kq — E = —y. Other 
interesting points become apparent, for instance, the one 
at y = — E = While there is otherwise no appar- 

ent relation to solutions of k s = +y, it should be noted 
that its s — case does produce a distinct feature in the 
line E = y. 



D. Case HN5 



In close correspondence with 



the treatment in 
Sec. IV B we can employ the RG in Eqs. (22) [23) [24) for 



transmission through HN5 consisting of a ring of N = 2 k 
sites. The sole difference with Sec. IIVBI is that all links 
of type Aj are initially present in this network. Yet, the 
details of the RG calculation in Sec. Mil show that under 
renormalization only links of type Ai renormalize while 
those Ai for i > 2 remain unrcnormalized at any step. 
The diagonal elements are again hierarchy-independent, 

k\ = E, while the recursion for A changes. Abbreviat- 
ing k = Ki, t = To, and A = Ai, 
their initial conditions reduce to 



Eqs. (22 23 24 1 and 



K m+1 
Am+1 



2A r , 

Am 
1 - 



(«o = E), 

(r = -l), (52) 

(A = -l). 



We have evolved the RG- recursion in Eq. ( 52 1 and plot- 
ted K ( fc = 10 ) as a function of k (0) = E in FigT^J At that 
(definitely not asymptotic) system size, N = 2 10 ss 10 3 , 
domains of localized states remain which will disappear 
asymptotically, as for the case of HN2 in Fig. [5[ Unlike 
for HN3, there are no steady-state solutions for Eq. (52) 
that could signal localization. It is interesting to ana- 
lyze the cause of this behavior. Since the A-links [green- 
shaded in Fig. [2] are an original feature of the network, 
the RG recursion in Eq. ( 52 1 for A m obtain a constant off- 
set preventing A m , and hence r m , from vanishing. Unlike 
for HN3, these small-world links allow perpetual trans- 
mission within a given level of the hierarchy, instead of in- 
terfering with other paths, and transmission is enhanced. 
In fact, we have studied also an interpolation between 
HN3 and HN5 by attaching a relative weight y to these 
A-links in HN5, relative to the otherwise uniform links 
present in HN3. Thus, for y = HN3 is obtained, and 
for y = 1 we get HN5. Yet, we find that for any y > 0, 
this model eventually behaves like HN5, with unfettered 
transmission throughout the energy spectrum. 



V. SUMMARY AND DISCUSSION 

We have devised a decimation RG procedure within the 
matrix methodology of Ref. [16] to obtain transmission 
of quantum electrons through networks within the tight- 
binding model approximation. This decimation RG pro- 
cedure of Appendix [X] can in principle be implemented 
for any network, in that Eq. ( Al ) with n + m sites is re- 



duced to Eq. (A7| with n sites. For general networks the 
bookkeeping required could become prohibitive. How- 
ever, we find the RG equations well suited to the networks 
we have chosen to analyze, namely the hierarchical net- 
works called HN3 (Hanoi Network with 3 bonds per site) 
and HN5 (Hanoi Network with an average of 5 bonds per 
site). This is because the RG can proceed block- wise, 
as depicted in Fig. [4[ We have analyzed both the ring 
geometry (Fig. [T) and the linear geometry (Figs, [I] and 
[3) of these networks, with the only difference being the 
last steps of the RG (Appendix B). We have also ana- 
lyzed how the transmission for the linear lattice (which 
we label HN2) changes with the strength of the small- 
world-type bonds added to form HN3 (Fig. [8) . 

The Hanoi networks are hierarchical models that pro- 
vide an intermediary between regular lattices and lat- 
tices that have random small- world bonds placed on reg- 
ular lattices. Since the small-world-type bonds in the 
Hanoi networks provide short-cuts between sites, one 
might expect intuitively that they should provide extra 
paths for transmission. However, because of the hierar- 
chical nature of the networks, the networks no longer 
possess translational symmetry. This broken symme- 
try is seen by the incoming quantum electrons, and can 
lead to Anderson-type localization. For the HN3 network 
with variable strength y for the small-world-type bonds, 
we find that the more we weight the small-world links, 
which classically would expedite transport [39], the less 
quantum transport is possible (Fig. [8) . Furthermore, we 
find that HN3 has band edges at particular energies E 
of the incoming electrons, between band gaps with near 
zero transmission and regions of extended wavefunctions 
and transmission near unity. The network HN5 adds still 
more small-world type bonds to HN3, but we find that 
for any non-zero strength of these additional bonds the 
band edges seen in HN3 disappear and approximately 
unfettered transmission is seen for large enough lattices 
for any energy of the incoming electrons. Thus the hi- 
erarchical nature of these lattices lead to very interest- 
ing transmission properties. Thus for these hierarchical 
lattices, the metal-insulator transitions depend on quan- 
tities other than just the embedding dimension. Similar 
effects have been seen for critical phenomena in hierar- 
chical lattices, but only where translational symmetry is 
broken [HI • 

Since the HN3 and HN5 networks are planar (see for 
example Fig. [3) , experimental realizations of these net- 
works should be possible to construct using etching tech- 
niques. These experiments would lead to very interest- 
ing device physics, in particular at the energy-dependent 
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FIG. 8: (Color online) Plot of the state of Eq. ( |50[ ) after 1000 iterations for initial energies — E and the interpolation 
parameter y (the resolution is 0.001 in each direction). Shaded points have not, or not yet, converged to a steady-state, i.e. 
those values possess non-zero transmission. At y = 0, corresponding to HN2, the system transmits for any input energy, 
— 2 < E < 2. As soon as y > 0, bands of localized states emerge (especially at E = 0), and the remaining transmitting states 
exhibit a chaotic dependence on the parameters. At y = 1, the case of HN3 marked by a dotted horizontal line, only a few 
non-localized states remain, and the further strengthening of small-world links diminish transmission even more, such that it 
ceases completely for y > 4. Any band of transmitting states appears to be accompanied by solutions of Eq. (51 1, which here 
are draw as curves for the lowest orders of the recursion only; a = 0, 1,2, and 3 corresponding respectively to dot-dot-dashed, 
dot-dashed, dashed, and solid curves. 
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FIG. 9: (Color online) Plot of k (10) in Eqs. g for HN5 as a 
function of its initial condition K*- ' = E. Bands of localized 
and delocalized states intermix, but those localized intervals 
disappear asymptotically. 



band edges we have analyzed for HN3. 



ing electron by E — 2cos(g). The reflection probability 



is R 



and the transmission probability is T = \t\ 



With this ansatz, the required solution of the infinite 
matrix Schodinger equation reduces to the solution of a 
finite matrix equation of dimension two larger than then 
number of sites in the 'blob'. Unlike Ref. [IB] , we assume 
no direct hopping between the leads, i.e. no short-cut 
path around the 'blob'. 

Consider the case with n + m sites in the blob. We 
specialize to the case where all hopping parameters (r or 
A) and on-site energies (convoluted with E to give k) are 
real. The (n + rn + 2) x (n + m + 2) matrix to solve for 
the transmission T = \t\ 2 is [IB] 



w 



d 



w 

'w A B it 
w d B T D u d 



u 



t 

Ipd 

V t J 



On 

V o J 



(Al) 



with the imaginary part of the complex function £, 
and the definition 
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Appendix A: General Matrix Formulation of the RG 



As in Ref. [TBJ, a 'blob' of atoms in the tight-binding 
approximation is considered to be connected to two semi- 
infinite leads. Each semi-infinite lead is considered to 
be a linear arrangement of tight-binding sites with on- 
site energy and a hopping parameter of —1 (setting 
the energy scale). The incoming (outgoing) lead is con- 
nected to the 'blob' sites by a vector of hopping param- 
eters W, (respectively, u). The Schrodinger equation for 
the infinite system, Wco^ = E^ must be solved with 
the appropriate boundary conditions. The ansatz made 
|16j is that the wavefunction at every site in the in- 
coming lead has the form tp m -i — e lmq + re"™ 9 with 
m = — oo, • • • , —2,-1, 0, and the outgoing lead has the 
wavefunction f/'m+i = te lmq with m = 0, 1,2, ••• , oo. 
The wavevector q is related to the energy of the incom- 



E 



E y/A-E 2 

h i 1 . 

2 2 



(A2) 



The matrix A is of size nxn and includes all interactions 
between the n sites that will remain after the RG. The 
matrix D is of size m x m and includes all interactions 
between the m sites that will be decimated by the RG. 
The matrix B is of size nxm and includes all interactions 
between the the n sites that will remain and the m sites 
that will be decimated. The matrices A and D are both 
symmetric matrices, while in general B is not symmetric. 
The vectors w, u, if), and n are all of length n, while the 
vectors Wd, u*d, tpd, and m are all of length m. 

Multiplying out the two middle rows gives the equa- 
tions 



and 



(1 + r)w + Aip + Bip d + tu = r , 



(1 + r)w d + B T 4> + T>ifj d + tu d = r , 



Solve Eq. (A4) for tp d to give 



^ d = -D 1 {l + r)w d + B T ^+tu d 



(A3) 



(A4) 



(A5) 



Substituting ip d into Eq. ( A3 ) and collecting terms allows the equation to be rewritten as 



[w — BD" 1 ^] (l + r)+ [A — BD _1 B T ] ip + [u — BD" 1 ^] t = 0„ 



(A6) 



Note that since D is symmetric, so is D . We can also substitute ip d from Eq. (A3) in for expressions obtained from 
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Linear geometry 



Ring geometry 
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Tl 




A| A, 



m = k — 2 




4 



/ 

i 
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«1 
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«1 r 



4 

TO 



«2 



'(1 



m = k — 1 



m = k 



w 



Ki 



2r n «i 



FIG. 10: (Color online) The last steps in the RG for the linear (left) and the ring (right) geometries. In both cases one is left 
with only two sites connected to the external leads. Note for clarity we have dropped the superscripts that denote the RG 
number on all variables. 



multiplying out the top and bottom rows of Eq. ( Al ). This gives that the matrix equation 



f i-wJ~D 1 w d V 
u;-BD ^ 
-uJjD- 1 



wjD _1 B T 





(A7) 



has the same solutions for r, t, and tp as does Eq. (Al). For the cases in this paper we will not have interactions 
between the input site and output site and the m sites to be decimated, so both Wd and Ud will be zero and the 
(n + 2) x (n + 2) matrix equation to solve for t becomes 



(A8) 



W 




BD 


^B T u 


II 







This completes the decimation RG of the m sites. 



Appendix B: Transmission from Small Renormalized RG steps require use of these equations with care. 
Lattices 



Equations! 20 - 24 1 can be used to construct the m — 
k — 2 state shown at the top of Fig. 10 The subsequent 
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1. Linear geometry 



For the linear geometry (left side of Fig. 10), for the 
to = k — 1 RG step we proceed as follows. Using the RHS 
of Eq. ( 20 1 with i — 2 we calculate 



~(fc-l) _ (fc-l) 

K 2 — K 1 



+ ftfc+i — Kfc — 2(A 1 



(fe-i) 



A fc ) (Bl) 



and then substitute n 



(fe-i) 



by k 



r,(k-i) 



in the RHS of 



Eq. (21) to get 



.(fc-i) 



i \ \( fe -!) 
K k+1 +\ k -\\ ' + 



(B2) 



This completes the to = k — 1 RG step. Now for to = k, we proceed in a similar way. Using the RHS of Eq. (22) 
with m = k — 1 and r k = 0, we calculate 



-(*) _ ffc-i) 



(fc-i) 



A fe ) 



2 k, 



(fe-l)l2 
3 J_ 

.(fc-l) 



(B3) 



and then with to = fc, i = 2, A^+i = and substituting k\ by in the RHS of Eq. (20), we calculate 



=(*) - =(*) 



K fe+2 — K fc+i — 2A^ ■* . 



(B4) 



Eq. (21), we get 



(k) (k) (k) 

Substituting to = k — 1 in Eq. (24) gives = 0. Finally, with m = k and substituting by k 2 in the RHS of 



Jfc) 



Kfe+i + A*; — A : 



(fc-i) [ r o 



(fc-l)l2 



.(fc-l) 



.(fc-l) 



= n k+1 + X k + r (fe) - 2Ai' 



.(fc-l) . !r (fc-l) 



(B5) 



After performing fc RG steps, we are left with just 1 + 2 k k = 2 sites with an on-site energy corresponding to re^ 

ffc) 

and an interaction of Tq ; between them. One site is connected to the input by w while the other is connected to the 
output by u. In order to decimate these two sites, we have 



o a ; 



D 



(fc) (fc)\ 
n Tc i . 
.(fc) ,» ' 



B 



A! = A BD B = 



w 
u 

(fc) 



(fc) 



(B6) 
(B7) 



where /3 = um/([/Ci ] 2 — [tq^] 2 ) and 7 = u/w. Thus after decimating all N = 1 + 2 fe sites, we have 



.(fc) 



(fc) 
(fc), 



1 + r 



2»9f(0 




From the above matrix equation by taking the inverse of the 2x2 matrix 



1 + r 



2*3(0 



(C-«i fe) /?/7)(^-<^7)-/? 2 [rn 2 



(fc), 



.(fc)i 



(B8) 



(B9) 



Thus we have found the transmission T = \t\ 2 . 



Now using Eq. ( A5 1 we can calculate the wavefunctions 



associated with the various sites. For simplicity we take 
w = u. To find the wavefunction associated with the 
highest level sites (i = k + 1), we start with the m = k 
RG step given in Fig. 10 Here, for Eq. (A5), ijj = 0, 



w d 



whereas the D and B matrices 



are given by Eq. Making these substitutions, we get 



ipk+i = 



(l + r)4 
(fc) 



.(fc) 



tr, 



(fc) N 
.(fc) 



(B10) 
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w d 



Here, for Eq. 

(fe-i) 
= K i 



(A5I, 



Next we can find the wavefunction associated with the 
immediate lower level site (i 

tp = V*H 
Sk-i) (1^ 



k). 
D 



= u d , D = kT l ' and B 

Making these substitutions, we get 



T (k- 
r o 



t) 



VI ft,-. 



-O) 



We can continue in this manner, in principle, to find the 
wavefunction of all the levels below. However the calcu- 
lations get very tedious and the expressions very compli- 
cated below i — k and therefore are not given explicitly. 



sites. Site 1 is the nearest neighbor of site and site 
2. This makes the interaction between the two sites to 



be 2r, 



(fe-i) 
o 



This adds a factor of 2A^ ^ to the on-site 



energy of site 0. Thus the on-site energy of the even site 
corresponds to (/tfe-i — 2lk-i) + 2Zfc_i = Kfc_i. The 
decimation of these two sites is similar to the case of 



the linear geometry, except that we need to replace k\ 



(k) 



(Bll) by 4 



(fc-i) 



and 

4[T (fe ~ 1) ] 2 ) and 



Ak) 



by 2t, 



(fc-i) 



So 



(fc-l)l 



2. Ring geometry 



Next consider the ring geometry (right side of Fig. 10 1 . 
We perform k—1 RG steps to be left with 2 k ~ ( - k -^ = 2 



1 + r\ = 2i^{0 / £ - ftf^/V 



(B12) 



This completes the RG, giving the transmission T = \t\ 2 . 

To find the wave functions here, we start at the m = k — 1 RG step and proceed as we did in the case of the linear 
geometry. We find that 



^(l + r ) K f- 1 )-2tT (fc - 1)N 
[K (*-i) ]2 _ 4[T (*-i) ]2 ^«C*- 1 )-2(l + r)r ( * 



and 



r< fc - 2) W (l + r + Q 



Notice that the two i = k—1 sites are symmetric. There- pressions for the lower i sites become complicated upon 
fore the wavefunctions obtained above for the two sites further iteration of this methodology, 
are identical, as expected. Again the wavefunction ex- 
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